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Estimating the Strength of an Elastic Network Using Linear Response 

Gemunu H. Gunaratne 
Department of Physics, University of Houston, Houston, TX 77204 
The Institute of Fundamental Studies, Kandy 20000, Sri Lanka 

Disordered networks of fragile elastic el ements have been p r oposed as a model for inner porous 
regions of large bones [Gunaratne et.al., :ond-mat / 0009221] , http://xyz.lanl.gov]. In numerical 
studies, weakening of such networks is seen to be accompanied by reductions in the fraction of load 
carrying bonds. This observation is used to show that the ratio F of linear responses of networks to 
DC and AC driving can be used as a surrogate for their strength. The possibility of using T as a 
non-invasive diagnostic of osteoporotic bone is discussed. 



Osteoporosis is a major socio-economic problem in an 
aging population Non-invasive diagnostic tools to 
determine the need for therapeutic intervention are es- 
sential for effective management of the disease. Bone 
Mineral Density (BMD), or the effective bone density 
is the principal such investigative tool [g]. Ultrasound 
transmission through bone H and geometrical character- 
istics of the inner porous region or trabecular architecture 
(TA) ji| are being studied as complementary diagnostics. 
In this Letter, we use results from a model system to sug- 
gest an additional diagnostic for osteoporosis. 

The TA is the principal load carrier in bone of older 
adults and cellular models have been used to study 
TAs H . In Ref. || , a disordered network of fragile elastic 
elements [Q was proposed as a system to model mechani- 
cal properties of a TA. Preliminary studies are conducted 
on two dimensional square networks which include elas- 
tic and bond-bending forces. Motivated by conclusions 
from mechanical studies of bone ||, the springs are as- 
sumed to satisfy a strain-based fracture criterion; specif- 
ically the fracture strains of elastic elements are chosen 
from a Wcibull distribution || (with parameters 7 e and 
to) . Bonds are assumed to fracture when changed beyond 
a critical angle; these fracture-angles are chosen from a 
second Weibull distribution (with parameters 7& and to) . 
Osteoporosis is modeled by random removal of a fraction 
v of springs from the network. 

The characteristic introduced below includes response 
of a network to an AC strain, which depends on the mass 
distribution on the network. It is modeled by placing 
masses m at the vertices. Viscous effects of the surround- 
ing medium are modeled by a dissipative force propor- 
tional to the speed of each mass. In studies reported here, 
points located on the sides were constrained to move ver- 
tically. Numerical studies of the system support the con- 
jecture that these elastic networks are a suitable model 
to study mechanical properties of bone H . 



Figure [l] shows the stress distributions on two net- 
works subjected to uniform compression 1 10 ll[] . For clar- 
ity only the compressed bonds are shown, and darker 
hues represent larger stresses. On the "healthy" network 
(where v — 10%) the large stresses supporting propaga- 
tion are seen to be distributed evenly over the network. 
In contrast, elastic elements supporting a "weak" net- 
work (y — 30%) form a few coherent pathways. We refer 
to the set active in load transmission as the "stress back- 
bone" of a network Jl2| . 




FIG. 1. The stress distributions on networks of size 40 x 100 
with (a) v = 10%, and (b) v = 30% representing "healthy" 
and "osteoporotic" bone respectively. For clarity only the 
compressed bonds are shown, and darker hues represent larger 
stresses. The crosses denote locations of long horizontal frac- 
tures. 

In Figure [j], the X's denote long horizontal fractures; 
specifically, locations where four or more consecutive ver- 



1 



tical bonds are missing. It is clear that these fractures 
prevent the participation of many bonds in the stress 
backbone. Since the number of long fractures increases 
with is, a progressively smaller fraction of bonds are able 
to be load carriers. The assumption that this dilution of 
the stress backbone is related to the ultimate (or break- 
ing) stress of the network motivates the measure intro- 
duced below. 

Consider first an ordered N x M square network of 
identical springs, each of whose breaking stress, break- 
ing strain and elastic modulus are denoted by Uo, Co and 
Y = Mo /Co- Then the ultimate stress and strain of the 
pure (i.e., v = 0) network are U(0) = Nuo and MCo- De- 
note by U {v) and C{v) their values for a network obtained 
by removing bonds with a probability v. Each horizontal 
layer of the network is approximately compressed by 
and a fracture will propagate when one of the bonds bor- 
dering it is stressed beyond uq. In this approximation, 
the DC response of the network is 
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(1) 



where u(v) = U(v)/U(0). 



in a network. First subject the network to a DC com- 
pression (Cdc) below the yield point, so that there is 
no fracture of elastic elements. Next introduce an ad- 
ditional AC compression, given by £(t) = CAC ex p(*^), 
where (ac "C Cdc @ ■ When f2 increases so does the at- 
tenuation of the signal, and a progressively thinner slice 
of the network is effected by the AC signal. Figure || 
shows the AC response of the network of Fig. |](b) driven 
at two frequencies. 

Denote by T(t) the sum of vertical forces on the top 
layer due to the AC strain. The linear response of the 
network x(f2) is given by T(Q) = x(Q) -C(^), where C(f2) 
and T(f2) are the Fourier transforms of C(i) and T(t) re- 
spectively. Figure || shows the behavior of |x(fi)| for a 
disordered network with v — 30%. When O — i> 0, there 
is no attenuation and hence x(fi) approaches xo- O n the 
other hand, for sufficiently large only the two edges are 
excited and hence x(£l) can be expected to approach a 
limit. Then each of the (1 — v)N bonds of the top layer 
are strained by a same amount and hence 



X (Q) w N(l - v)Y, 



(2) 



where Y denotes the mean elastic constant for bonds on 
the top layer. Numerical integrations of disordered net- 
works confirm these conclusions, see Fig. ^. 
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FIG. 2. The distribution of stresses on the disordered net- 
work of Figure 1(b) due to small amplitude AC compressions 
with (a) Q — 10 and (b) Q = 500. As Q increases, so does the 
attenuation of the signal, and a progressively thinner slice of 
the network is effected. For sufficiently large Q only bonds 
belonging to the top and bottom layers experience an AC 
stress. 

Next we argue that the linear response to an external 
AC strain can be used to estimate the number of bonds 




FIG. 3. The dependence of |x(^)l ° n ^ f° r a network with 
v — 30%. For r2 — > 0, — > xo, while for sufficiently large 
57, x(fi) ~ N(l — v)Y. In this limit, only bonds bordering the 
top and bottom layers experience AC stress. The fluctuations 
for intermediate values of il are configuration dependent. 

Consider once again the ordered pure network. When 
only the top and bottom layers experience AC stress the 
effective height of the layer is reduced by a factor M; 
thus, x(fl) will be larger than xo by a factor M. For 
v = 0, this can be expected to hold (approximately) 
even for disordered networks since the stress backbone 
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covers the entire network. Even though both xo an d 
x(fi) will decrease with increasing v 1 the latter will be 
effected less because long fractures in the middle of the 
network (which reduce the stress backbone) have no ef- 
fect on x(O). These arguments motivate the use of 



I» = lim M- — 7~tt 



(3) 



to estimate the reduction in the extent of the stress back- 
bone. Approximating Y by Y and using Eqs. ([!]), (g) and 
© gives 



«(!/)« (i-i/).r(i/). 



(4) 



To complete the derivation, we propose a relationship 
v = F(u). Since u(u) for a 2D square network vanishes 
at the percolation threshold v = uq = h, 



= F(0). 



(5) 



A form for Fix) has been presented in Ref. JJ3J], where 
it was assumed that the propagation of the longest hor- 
izontal fracture causes the collapse of a network. The 
ingredients used are (1) a horizontal fracture of size k en- 
hances the stress on the bordering elements by (1 + ak a ), 
and (2) the size k m of the longest fracture satisfies 
MNv km {\ — v) 2 |l4| . Since the possible two dimensional- 
ity of a fracture is ignored, the resulting expression does 
not satisfy Eq. (||). We propose a modification 
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which agrees with results from numerical integration of 
disordered networks for v G (0,0.4) Jl5| . [T6[ | . Parameters 
a and a are expected to depend on factors such as the 
relative strengths of the clastic and bond-bending forces. 
Using Eqs. (@) and §), 



u{v) = (i-^ )r + /i(r), 



(7) 



where the nonlinear correction h(T) can be obtained by 
inverting (e.g., Born expansion) Eq. (^) and will depend 
on a, a, and ln(MiV). 

To test the validity of Eq. (fj]), numerical studies were 
conducted on a group of five equivalent disordered net- 



works 1 10 1 . T(v) and U{v) for a given network is evalu- 
ated using methods discussed earlier pHp^| . As can be 
expected, the normalization of U{v) by {7(0) reduces the 
variability between distinct configurations in the group. 
In evaluating the linear response, the signal T(t) was 
collected after the transients have settled. For a given 
v, there is scatter in the values of both T(v) and u(v), 



and Figure ^ shows their mean and standard errors. The 
dashed line is the best fit to the data in the form (Q) 
with h(T) — cT@. Numerical studies on configurations 
with different parameters show that, unlike the linear 
term of Eq. (ffl), the coefficients (3 and c are parameter 
dependent. 
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FIG. 4. The relationship between the function T(v) and 
u(v) for a set of 5 equivalent disordered networks. The dashed 
line shows the best fit of the data to Eq. (0) with h(F) = cF . 
For these parameters [10], c = 0.90 and (3 — 3.20. 

The reduction of bone strength from its peak value de- 
termines the level of osteoporosis. Unfortunately, it is 
not accessible in- vivo (without breaking a bone!), and 
surrogates such as bone density are used to identify os- 
teoporotic bone. The BMD of a patient is compared 
with that of a sample population to determine if and 
when therapeutic interventions are necessary. However 
the ultimate stress is known to depend on other factors 
of bone including the structure of its TA and the "qual- 
ity" of bone material. The resulting variations make it 
difficult to identify individuals susceptible to fracture us- 



ing measurements of BMD alone 17 1. 

In this Letter, we have shown that linear response of a 
network can be used as a surrogate for its strength. This 
conjecture is based on the assumption that the strength 
is related to the extent of the stress backbone. Mea- 
surements required to evaluate T can be implemented 
on ex-vivo bone samples. DC strain can be imposed us- 
ing pressure loading, and protocols using ultrasonic tech- 
niques have been developed to evaluate the response of 
bone samples to AC driving [Q. Previous studies sug- 
gest that when their frequencies larger than ~ l.bMHz 1 
ultrasonic signals will excite only those trabeculae on the 
outer edges of a TA [^9|. How these measurements can 
be implemented in-vivo remains to be studied. 
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Several issues need to be reiterated. To calculate x(O), 
only elastic forces (on vertices of the top layer) were 
used. In driving a bone sample with an AC strain, the 
matter (on the outer layer) is accelerated in a dissipa- 
tive medium. These inertial and dissipative forces are 
proportional to Q 2 and fl respectively. In contrast, for 
sufficiently large fl, x(O) is Sl-independent (see Fig. ||). 
Hence, the latter can be extracted from the response of 
the TA to AC signals of several frequencies. Secondly, 
observe that for weak networks (y close to Uq) the re- 
lationship between u(v) and T(is) depends only on vq. 
For smaller values of v, however, the nonlinear correction 
h(T), which depends on model parameters, becomes rele- 
vant. Hence the form of u(T) may have to be determined 
for distinct bone locations (e.g., femur, vertebrae) before 
a complete diagnostic tool for osteoporosis is developed. 
Finally, notice that the definition of T includes the num- 
ber of layers M of the network. Since the lengths of 
trabeculae from specific anatomical locations are known 
(typically ~ 1mm) , the length of the sample can be used 
to determine M. 
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